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Abstract. Recently, the issue of the role of the Eddington limit in accretion discs became a mat- 
ter of debate. While the classical (spherical) Eddington limit is certainly an over-simplification, 
it is not really clear how to treat it in a flattened structure like an accretion disc. We calculate 
the critical accretion rates and resulting disc luminosities for various disc models corresponding 
to the classical Eddington limit by equating the attractive and repulsive forces locally. We also 
discuss the observational appearance of such highly accreting systems by analyzing their spectral 
energy distributions. Our calculations indicate that the allowed mass accretion rates differ con- 
siderably from what one expects by applying the Eddington limit in its classical form, while the 
luminosities only weakly exceed their classical equivalent. Depending on the orientation of the 
disc relative to the observer, mild relativistic beaming turns out to have an important influence 
on the disc spectra. Thus, possible super-Eddington accretion, combined with mild relativistic 
beaming, supports the idea that ultraluminous X-ray sources host stellar mass black holes and 
accounts partially for the observed high temperatures of these objects. 
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1. Introduction 

With the constant improvement of observational techniques in the past decades, more 
and more detailed information about accretion disc systems could be gained. In par- 
ticular, a large number of ultraluminous X-ray sources (ULXs) have been discovered 
(Fabbiano 1989), imposing a severe problem upon the existing general idea of accretion 
disc systems: With a bolometric luminosity exceeding lO'^^ergs"^, at least some of them 
show relatively low radiation temperatures 0.1 keV). These systems have been sug- 
gested to be intermediate mass black hole (IMBH), sub-Eddington accretion disc systems 
(Roberts, Warwick, Ward, et al. (2005)). However, large samples of ULXs from recent 
observational data reveal that a distinct class exists, showing higher temperatures - some- 
times exceeding 1 keV - than can be explained by IMBHs (Mizuno, Ohnishi, Kubota, 
et al. (1999)). In contrast, stellar mass black holes accreting above their Eddington limit 
can account for these sources (Watarai, Mizuno & Mineshige (2001)). 

In particular, the applicability of the classical Eddington limit is an important point 
in this discussion: Similar to the stellar case, disc accretion may be limited by radiation 
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pressure, counteracting gravity and viscous dissipation. While in the stellar case, we 
are dealing with an approximately spherically symmetric, i.e., 1-dimensional situation, 
discs require an - at least - 2-dimcnsional treatment. The stellar Eddington limit implies 
several assumptions, which do not apply for the disc case, though it is not clear a priori 
to what degree a proper treatment will alter the resulting numbers. 

In our work, we calculate the local analog to the Eddington limit for a classical thin 
disc model of the Shakura & Sunyaev (1973) variety and for a slim disc model of the 
Abramowicz, Czerny, Lasota, et al. (1988) type (Sect. 2). In both descriptions, a standard 
a-viscosity is applied, therefore restricting the models to the non-selfgravitating case. 
We also investigate the spectral energy distribution of supercritical accretion flows in the 
stellar mass black hole case, based on the radiation hydrodynamic (RHD) simulations 
computed by Ohsuga, Mori, Nakamoto, et al. (2005) (Sect. 3). 

2. Black hole accretion and the Eddington limit 

The physical reasoning for the Eddington limit is an equilibrium of the attractive 
force Fg and the repelling force F^. The original stellar Eddington limit relies on several 
assumptions: (1) spherical symmetry of the system; (2) isotropic radiation; (3) homo- 
geneous degree of ionisation; (4) Thomson scattering as the sole source of opacity; (5) 
negligible gas pressure; (6) no relativistic effects; and (7) no time dependence (station- 
arity). To investigate if these assumptions are an oversimplification in the disc case, we 
relinquish the approximations (l)-(5), but keep the approximation of non-relativistic 
stationarity for our calculation. 

We assume azimuthal symmetry in the following. Therefore, one can define two differ- 
ent "Eddington limits" , corresponding to the equilibria in the vertical {z) and radial (s) 
direction.'^ For a rotating viscous disc, further contributions to the total force have to be 
included in the radial direction: 

F^:l {s) = F W is) - F W {s)=0, F W (s) = F W (5) - F^^^ (s) + . . . = . 

In the following, as an example, we calculate the maximum amount of matter Merit 
that can be accreted towards the central object in the case of a stellar mass black hole, 
M = 10 M0, for non-selfgravitating discs with an a-parameter of 0.1. For both the thin 
and the slim disc case, we use two descriptions for the opacity: Firstly, simple Thomson 
scattering is applied: 

'^es — 0*4 cm g . Secondly, we use an interpolation formula, k — Kin 
(Gail, priv. comm.), which accounts for further contributions like bound-free and free-free 
absorption. For a detailed description of Heinzcller & Duschl (2006). 

We refer our results to the classical Eddington- limit, which in the disc case with a 
torque- free boundary condition at the disc's inner radius Sj = 3rg = 6GM/ is given by 
Le = 4ncGM/Kes = 1.2 ■ lO^^erg/s, and Me = 87tcsi/Kes = 2.6 • 10"'^ Mo/a. Detailed 
calculations reveal that the vertical limit sets severer conditions on the critical accretion 
rates and corresponding luminosities than the radial limit in both the thin and the slim 
disc model. We therefore concentrate on the results for the vertical limit in the following. 

Figure 1 shows the results for the critical accretion rates Merit in the vertical Eddington 
limit as a function of the radial distance s from the central object. In the slim disc case, 
one free parameter remains to be determined by hand: the ratio of the local height h of 
the disc to the local radial distance s. Since the slim disc model is valid for h/s < 1, we 
plot the results for the two cases (/i/s)max = 1-0 and (/i/s)niax = 0.5. 

Obviously, the critical accretion rates are no longer given by a global quantity like Me- 

^ We use a cylindrical coordinate system {s, ip, z} with the distance r to the origin, = s'^+z'^ 
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Figure 1. Critical local accretion rates in 
the vertical Eddington limit as a function 
of the radial distance s from the central ob- 
ject. For slim discs, (/i/s)max = 1-0 in the 
upper and (/i/s)max ~ 0.5 in the lower case. 
Black lines correspond to Thomson scatter- 
ing, K = Kos, gray lines to interpolated opac- 
ities, K — ^iii 
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Figure 2. (a) Spectral energy distribution for various inclination angles B; (b) Comparison of 
total intensity, average photon energy and total photon number for = 0... 7t/2, normalized 
to the corresponding values at 7t/2 



they depend on the radms oc si-2- -i-9^ While the opacity has a strong influence on the 
results (A/crit oc k^^) for low accretion rates, it becomes more and more unimportant 
in the slim disc model for high accretion rates. This is due to the fact that advection 
dominates the energy transport in this regime. As can be seen from the results for thin 
discs, the inner boundary condition plays a crucial role in determining the maximum 
amount of matter which finally reaches the central black hole. For slim discs, the critical 
accretion rates may reach up to IO^Me, while at the same time the luminosity stays close 
to the classical limit, L « 20Le (Watarai, Fukue, Takeuchi, et al. (2000)). Thus, highly 
super-Eddington accretion comes along with only mildly super-Eddington luminosities. 



3. Spectral energy distribution of super-Eddington flows 

In the second part, we investigate the observational appearance of super-Eddington 
accretion discs. This work is based on 2D RHD simulations, which are presented in detail 
in Ohsuga, Mori, Nakamoto, et al. (2005). In a spherical computational box of 500i?G 
radius, matter is injected continuously in the disc plane x — y around a IOMq black hole 
at high rates [M — IO^A^e)- Again, the viscosity is parameterized by a = 0.1. For the 
energy transport equation, radiation and advection are taken into account, as well as 
photon trapping effects (Begelman (1978)). We calculate the spectral energy distribution 
of this system for various inclination angles 8 between (face-on) and 7t/2 (edge-on). For 
details about the computation, we refer the reader to Heinzeller, Mineshige & Ohsuga 
(2006). The resulting SEDs for inclination angles 0, 7t/4 and 7t/2 are shown in Fig. 2(a). 



4 



HeinzcUcr, D., Duschl, W.J., Mincshigc, S. & Ohsuga, K. 



While in the Rayleigh- Jeans part of the spectrum no difference between the orientations 
of the system can be seen, the peak intensity region around 1 keV is shifted towards higher 
photon energies and higher luminosities for lower inclination angles. For all orientations, a 
plateau-like structure in the high-energy part of the spectrum can be observed, although 
this feature is more prominent for the face-on case. Its origin lies in the emission of the 
hot gas in the vicinity of the black hole: when calculating the spectrum for = while 
neglecting gas emission and absorption, this plateau disappears. 

To outline the orientation effects, we compare in Fig. 2(b) the total intensity emitted by 
the system for inclination angles 8 = {0 . . . 7t/2}. For small inclinations, the intensity rises 
by a factor of up to 6.4 for = 0; this is due to both an enhanced average photon energy 
and an increased total photon number. We identify this behavior with mild relativistic 
beaming. When fitting the peak intensity region with a black body spectrum, we derive 
temperatures between 1.1 • 10*^ K for 9 = 7t/2 and 1.7 • 10*^ K for 6 = 0. 

4. Conclusions 

The above investigations allow us to draw several conclusions on black hole accretion 
processes. Firstly, wc show that the classical Eddington limit can not be applied in 
discs. Contrary, the critical accretion rates become a local quantity, depending on the 
distance from the central object. Thereby, the inner disc region and inner boundary 
hold the key position in determining the final amount of matter that can be accreted 
by the central object. We find that super-Eddington accretion is possible with accretion 
rates up to IO^Me, corresponding to slightly super-Eddington luminosities up to 20Le- 
Since the bottle-neck for accretion lies in the inner region of the disc, our results may 
also correlate to the question on the origin of outflows and jets. Secondly, our work on 
the spectral energy distribution of such super-Eddington accreting systems shows that 
for a proper interpretation of the spectrum, the combined system of the disc and its 
surroundings has to be taken into account. Features like the high-energy plateau can not 
be reproduced by considering the accretion disc only. Also, we find that mild relativistic 
beaming becomes important for small inclination angles, leading to higher luminosities 
and also higher temperatures. However, these temperatures are still too low to account 
for the observed ULX sources which sometimes reach temperatures up to 10^ K. This 
may be due to the lack of Comptonization effects in our calculations: Socrates, Davis & 
Blaes (2004), for example, showed that turbulent Comptonization produces a significant 
contribution to the far-UV and X-ray emission of black hole accretion discs. Thus, the 
inclusion of Comptonization effects will be a primary goal in our future investigations. 
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